Imaging spin-inelastic Friedel oscillations emerging from magnetic impurities 
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We consider inelastic scattering of localized magnetic moments coupled with the electrons on the 
surface. We argue that spin-inelastic transitions of the magnetic impurities generate oscillations at 
a momentum k, corresponding to the inelastic mode, in the second derivative of the current with 
respect to voltage d 2 I/dV 2 . These oscillations are similar in nature to Friedel oscillations. Inelastic 
Friedel oscillations, which were previously proposed for spin-unpolarized set-up, is here extended for 
spin-polarized systems. We propose to use scanning tunneling microscope to measure spin-inelastic 
scattering generated at the impurity spin by imaging the d 2 I/dV 2 oscillations on the metal surface. 

PACS numbers: 73.63.Rt, 07.79.Cz, 72.25.Hg 



Measurements of inelastic transitions open up a route 
to investigate the excitation spectrum of physical sys- 
tems. There has been a growing activity in elucidating 
inelastic scattering processes in quantum systems using 
various experimental techniques. A n in complete list in- 
cludes inelastic neutrorP^l anc | X-ray^^ scattering, trans- 
port through break junctions^', and scanning tun neli ng 
microscopy (STM) with spin- polarized (SP-STM^H 

or 

non spin-polarized tipPHH. 

Surface imaging of scattering states can be performed 
e.g. by using STM to probe the spatial spectral den- 
sity variations at a given energy. It is well-known that 
Friedel oscillations emerge around defects ads orbed onto 
a surface caused by elastic scattering processe d 16 * 17 !. Less 
known is the prediction made by us earlier that points 
to the existence of inelastic Friedel oscillations emerging 
from vibrating impurities 1 These oscillations were re- 
cently demonstrated experimentally for dimers of meta- 
dichlorobenzene^. The mechanism for inelastic Friedel 
oscillations is essentially the same as for conventional os- 
cillations and comes from interference of incoming and 
outgoing waves that have an energy mismatch given by 
energy transferred to/from local vibrational mode. 

In this paper, we propose a spin-polarized extension of 
inelastic Friedel oscillations that arise from spin-inelastic 
transitions. We also propose to use STM to image these 
oscillations. A local magnetic moment interacting with 
surface electrons, generates a local spin-polarization in 
the surface. The spin-inelastic transitions provide an- 
other modification to the electronic and magnetic struc- 
tures of the surface states. This local modification ex- 
perimentally can be stimulated by adding or remov- 
ing an energy quantum that corresponds to the inelas- 
tic transition energy. Such an energy change causes a 
change in the electronic structure of the type 6N(r, w) = 
J2 a p 2a/3( r i a; )$( a; — Ap a ), for low temperatures. Here, 
Ap a is the transition energy, whereas Q a p(r,u>) is a spa- 
tial distribution function which depends on the involved 
states \a), \(3). 

In order to enhance the signature of spin-inelastic scat- 
tering effects, and image spin-inelastic Friedel oscilla- 



tions, we propose to use quantum corrals constructed out 
of magnetic atoms, or molecules. Placing STM tip at the 
center of quantum corrals will allow to amplify the signal. 
As the magnetic atoms are coupled through exchange in- 
teractions, the magnetic structure of the quantum corral 
can be engineered to meet specific requirements. When 
coupled anti-ferromagnetically, one can study qualitative 
differences in quantum corrals comprising an even or odd 
number of atoms. Ferromagnetic coupling between the 
atoms, on the other hand, gives rise to a large collec- 
tive spin moment that could have its own signatures. In 
either case of ferro- or antiferromagnetically coupled cor- 
rals using STM would allow to image specific spatial fin- 
gerprints throughout the interior of the corral. 

In order to give an example of the effect we are propos- 
ing, consider a collection of general (quantum) spins 
S„ = S(r n ) located at the positions r„ on a metallic 
substrate surface. The surface electron density can be 

modeled by "H S urf = J2w £ kcj ccr Ckcr, whereas the Kondo 
interaction between the surface electron density and the 
local spin is given as Hk — v u Jk J2 n s ( r «) ' ^n, where 
v u is the unit area and Jk is the Kondo exchange pa- 
rameter, whereas s(r„) = cj 7 (r n )cr . CT /c cr '(r„) with spin 
indices a, a' =t,|, c CT (r) = J c ko .e ikr dk/(27r) 2 , and the 
vector er of Pauli matrices. It should be noticed that 
the Kondo interaction only provides the isotropic inter- 
action between the local spin moment and the substrate 
electrons, whereas the anisotropy is being treated sepa- 
rately, see the discussion below. 

We make contact with current STM measurements on 
local magnetic moments by formulating the tunneling 
(differential) conductance in terms of the Terso ff and 
Hamann approach^, and its generalizations^^!. The 
tunneling conductance at low temperatures is in this ap- 
proach given by 

d ^Xl^n{e F -eV)N{v 1 e F ) (1) 

where Ef is the Fermi level of the system in equilibrium. 
This expression, thus, relates the tunneling conductance 
to the electronic (N) structure of the substrate surface, 



2 



and correspondingly (n) for the tip. It is, therefore, suf- 
ficient to study the local variations in the spin-polarized 
surface density of electronic states. 

In our calculated examples below, we assume that the 
localized spin moments can be described in terms of the 
Hamiltonian 

U S = Y^WSn? + E{(S+) 2 + (^) 2 ]/2}, (2) 



where the anisotropy fields D and E account for the effec- 
tive interaction between the localized spin moments and 
the surface electrons. Notice that this Hamiltonian ef- 
fectively describes the interactions between the localized 
spin moment and the surface electrons which give rise 
the anisotropy of the localized spin. The isotropic inter- 
actions are accounted for by the Kondo interation "Hk- 
This model defines 25" + 1 eigenstates |a) and eigenener- 
gies E a , and we introduce the operators d' a (d a ) which 
create (destroy) a particle in the state \a). For later use, 
we also define the spin operators rig = \a)(a\S l \/3)(/3\, 
i = x,y,z. 

We employ the model given in Eq. pi) since it has 
been successfully used to describe single (and multiple) 
impurities located on metallic surface, see e.g. Refs. ITOl - 
1141231 The anisotropy fields D and E are related to the 
properties of the interactions between the local adsor- 
bant and t he su bstrate material, and can be fitted to the 
exp eriment! 1 1 1 1 1 but also determined through first princi- 
ples calculations^ 

In absence of impurities, the substrate surface is as- 
sumed to be non-magnetic, however, in presence of the 
local spins the surface LDOS may become spin-polarized 
locally around the spins. We account for the scatter- 
ing off the magnetic impurities by calculating the real 
space Green function (GF) for the surface electrons us- 
ing G(r,r';zw) = J G(k, k'; iw)e lk r - ik ' r 'dkdk7(27r) 4 , 
where G(k, k'; iui) — \G aa i (k, k'; iui)}^ is the 2x2 ma- 
trix of the spinor *(k) = (ckf c^) T . 

First we construct a bare GF G^ ^ which contains the 
spin-polarization induced by the localized magnetic mo- 
ments, using the model Hsmf+H-K- In spirit of scattering 
theorj£4j25|, wg q^^^ a T-matrix (spin space) formula- 
tion 

G(°)(k,k') =<J(k - k')g(k) 

+ £ g(k) e - 4k - r "T(r„, r m )g(kV k '- r "* , 



(3a) 
(3b) 



where t _1 (r„,r m ) = <5(r„-r m )-V„g(r n -r m ), andg(r- 
r') = Jg(k)e ik -( r - r ')dk/(27r) 2 , g(k;iu;) = (ioj - e k )-\ 
The scattering potential V = Vq + <r • A n (iui) comprises 
the spin-independent and spin-dependent contributions 
V and A n (icj) = v u J K (S n )(iu). 

The spin-inelastic scattering off the localized spin mo- 
ments that influences the surface electron GF, is ac- 



counted for in second order perturbation theory with re- 
spect to v u Jk- We define the self-energy^ 



^trcr' ( r n j I'm i 



(VuJk) 2 



x G^](r„,r m ;zw - iv) 



(4) 



(/?-! = k B T) where the spin-spin GF Xnm( z ) = 
)'(-i)(TS rl (<)S m (t'))e lz(t - t ' ) dt', and obtain the real 
space GF 

G(r,r') wG (0) (r,r') 

+ G(0) (r, r„)S(r„, r m )G^ (r m , r'). (5) 



Next, we discuss the effects of the spin-inelastic scat- 
tering on the surface electrons. For the sake of argument 
we use the unperturbed surface electron GF, i.e. replace 
G(°) by g in Eq. to allow for analytical calculations. 
For non-interacting spin impurities, the spin-spin GF can 
be written 



T aj3 Ti3 a G n p{iu] + iv)G na (iv) (6) 



where G na (t,t') — (— z)(Td na (t)d^ a (f')). Here, we have 
written the spin in terms of the eigenstates of Eq. |2|) 
such that the spin GF G na (iuj) = (iu> — Ena)" 1 . We find 
that the retarded form of the self-energy in Eq. ^ can 
be written (using quadratic dispersion = h 2 k 2 [2m) 



a/3 

+ \f(E p ) - f(E a )] 
+ nrf(uj -E a + E f3 



In 



\u-E a +Ep\ 



W 



(7) 



where 2W ~ f eV is the band-width, whereas 7 = 
v u JxN a , with the bare surface DOS N — m/ti 2 , 
and f(x) is the Fermi function. Noting that g r (r) ps 
— iNo Jo(fe|r|)/2, where Jq(x) is the zeroth order Bcssel 
function of the first kind, whereas k = |k| = y^2Nouj, 
we find that the substrate LDOS N(r,u) — N n (r,uj) + 
SN(t,lj), defined by N(t,lj) = — tr ImG(r, r)/%, where 
JVo(r,w) = Nq whereas 

S -^p^=^J 2 (k\r-r \)J2T a p-T f3a {f(-E a )f(E p ) 

a/3 

+ [f{E^)-f{E a )]f{uj-E a + Ep)). (8) 

It is clear form this expression that the amplitude of the 
inelastic signal scales with the square of the Kondo cou- 
pling J K . 
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FIG. 1: (a) DOS (dl/dV) map of a single spin 5=1 impu- 
rity adsorbed onto a metallic surface around which concentric 
Friedel oscillations emerge, calculated using the GF defined in 
Eqs. and {H]). (b) IETS (d 2 I/dV 2 ) spectrum, with (solid) 
and without (dashed) spin-inelastic scattering, recorded at a 
distance 1.8 nm from the impurity (star in panel (a)). The 
inset shows the IETS spectrum at the defect, where the two 
plots correspond to the broadened (bold) and non-broadened 
(faint) spin states, (c), IETS maps of the system different en- 
ergies, marked by bullets in panel (b). Here, D — —10 meV, 
£~ \D\/5, T~4K. 



Feeding the energy to = E a — Ep into the system 
by means of e.g. the bias voltage, stimulates the in- 
elastic spin transition |a)(/3|, and the onset of the in- 
elastic scattering generates an abrupt change in the 
surface LDOS. The expression in Eq. ffify, moreover, 
shows that the onset of the spin-inelastic scattering gen- 
erates spatial variations in the charge density emerg- 
ing from the localized magnetic moment, referred to 
as spin-inelastic Friedel oscillations, in analogy to pre- 
viously introduced inelastic Freid el oscillations emerg- 
ing from vibrational defect s^ 1 ^. The charge density 
variations are modulated by the momentum k. Gener- 
ally, inelastic scattering is not a Fermi surface effect but 
rather pinned to the momentum fc, for which reason one 
should expect a varying wavelength of the Friedel os- 
cillations emerging from the impurity depending on the 
energy of the specific inelastic transition. The LDOS 
connects with current inelastic electron tunneling spec- 
troscopy (IETS) measurements by noting that the signal 
d 2 I(r, V)/dV 2 oc dN(r,uj) /du. For low temperatures, 
the derivative —df(u — E a + Ep)/duj -> <5(cj - E a + Ep), 
which indicates the possibility to image inelastic Friedel 
oscillations in a narrow range of energies around the in- 
elastic transfer energy E a — Ep. 

In Fig. [l](a), we plot the local DOS of the surface elec- 
trons interacting with a localized S = 1 spin moment ad- 
sorbed onto the surface, pertinent for e.g. Co/Pt(lll) 14 , 
around which elastic Friedel oscillations emerge in the 
surface DOS. For the calculations we have used the full 
electronic GFs prescription as defined in Eqs. ([3| and 
([5]). Positioning the STM tip at the point marked by 
a star in panel (a), we plot in Fig. [I] (b) the IETS 
(d u N(r t i P) u))) spectrum for the perturbed (solid) and 



unperturbed (dashed) surface. Here, we have added a 
phenomenological Lorentzian broadening (~ 7.5 meV) of 
the spin states in order to capture the behavior of the 
IETS spectrum observed in RefP^. The broadening has 
been estimated the by fitting the shape of the IETS spec- 
trum at the defect position to experiments^™ see inset of 
Fig. [l] (b) (bold). The faint plot in the inset of Fig. [I] 
shows the IETS spectrum for the spin in the atomic limit. 

The setup, thus, demonstrates the possibility to re- 
motely record the inelastic signatures emerging from the 
scattering center, due to its propagation over the surface 
via the spin-inelastic Friedel oscillations. 

The spatial characteristics is expected to vary signifi- 
cantly with the energy, which indeed can be seen in Fig. 
[I] (c), where we plot IETS maps for a few energies cor- 
responding to the energies in Fig. [I] (b) . As is indicated 
in the IETS spectrum, no essential spatial structure is 
found for energies far off the inelastic transition energies. 
In fact, since the IETS spectrum is vanishingly small for 
energies off the inelastic transition energies, the spatial 
IETS maps are expected to be nearly equal to the cor- 
responding bare maps. By a comparison between the 
IETS maps for the different energies, it is clear that the 
localized moment generates a spatial response, i.e. spin- 
inelastic Friedel oscillations, for probe energies close to 
the inelastic transition energies. 

The magnetic structure M(r,w) emerging from the lo- 
cal spin moment, is in this simplified example reduced 
to a simple spatially varying spin-polarization M(r, uS) = 
M z (r,u)z of the surface electrons, which can be calcu- 
lated from M z (r,w) = J2a a aa^<y{^i w). For weak cou- 
pling between the localized spin and the electron medium 
assumed here, the spin-polarization is negligible. 

We now consider quantum corrals comprised of mag- 
netic atoms and implement the theoretical framework in- 
troduced above. In particular, we focus on circular ge- 
ometry and consider the resulting electronic structure for 
independent spin moments. In Fig. [2] (a) and (d) we plot 
the Fermi level spectral density of circular and ellipti- 
cal quantum corrals comprising 20 independent S = 3/2 
atoms anti-ferromagnetically coupled. The structure in 
the spectral density is caused by the confinement imposed 
by the corral. The corresponding dl/dV calculated at 
the center of circular corral s with different radii, using a 
level broadening of ~ 5 me V 14 * 23 ^, is shown in Fig. [2] (b) . 
Due to the voltage dependent corral dl/dV, we expect 
d 2 I(r c , V)/dV 2 to be non-zero also for voltages far off 
the spin-inelastic modes, which is shown in Fig. [2](c). 

For the elliptic corral, in Fig. [2](e) we show the dl/dV 
at the focal point — c, c.f. labels in Fig. [2] (d), for cases 
with (bold) and without (faint) a spin defect at the focus 
c. It is clear that the electronic structure changes slightly 
due to the additional defect. More important is that 
the dl/dV becomes distorted near zero bias voltage due 
to the spin-inelastic scattering, and those distortions are 
more clearly seen in Fig. [2] (f), where the corresponding 
IETS spectra is plotted. Despite the energy variations of 
the surface electron density, which accordingly are also 
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FIG. 2: (a), (d) dl/dV map of a 20 atom circular and el- 
liptical quantum corral (atomic positions marked by pen- 
tagons), respectively, for independent S — 3/2 spins anti- 
ferromagnetically coupled, (b) dl /dV calculated at the center 
of circular corrals for radii 7? = 39.75 + {0.2, 1.2, 2.2} A, and 
(c) corresponding IETS (d 2 I/dV 2 ) spectrum. Inset shows 
the IETS for a single S = 3/2 defect, (e) and (f) dl/dV 
and d 2 I/dV 2 , respectively, of the corral in (d) calculated at 
(x, y) — (— c, 0) with (bold) and without (faint) an adatom at 
(x,y) = (c,0). Ellipse in (d) given by R 2 = (x/a) 2 + {y/bf, 
with R = 24.25 A, and a/b = 1.5. Here, D = -3.25 meV, 
E = 0, and T = 4 K. 



picked up in the IETS signal, the spin-inelastic contribu- 
tion provides a significant distortion of the signal. Our 
calculations performed for the quantum corrals, hence, 
provide clear demonstrations that the inelastic scatter- 
ing should be remotely detectable also within systems 
with more complicated electronic structures. 

The slow energy dependence of the underlying DOS 
implies that its energy derivative is small, which leads to 
that the spatial signatures in the IETS maps are corre- 
spondingly small for energies sufficiently far off the in- 
elastic transition energy. Close to the inelastic transi- 
tion energies, on the other hand, we expect to be able 
to detect the spatial variations in the spectral density, 
analogous to the maps shown in Fig. [T] (c). 

We have demonstrated theoretically, that it should be 
possible to image the response to spin-inelastic transi- 
tions using STM for IETS measurements. Scattering off 
the local spin (or magnetic) moment modifies the DOS 
locally around the impurity, and for specific energies cor- 
responding to the inelastic transition energies, additional 
modification of the local DOS is expected. The inelastic 
signatures can be identified as sharp peak/dip features 
near the inelastic transition energies. Performing IETS 
measurements will reveal a spatially modulated spectral 
density near the inelastic modes. While measurements 
on single magnetic impurities should be sufficient in or- 
der to resolve the inelastic Friedel oscillations, we suggest 
that the corresponding signatures can be substantially 
enhanced by engineered quantum structures, e.g. quan- 
tum corrals. 
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